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Abstract 

A general solution for the electrostatic potential in an atomic-thin-body (ATB) field-effect 
transistor geometry is presented. The effective electrostatic scaling length, Xeg, is extracted from 
the analytical model, which cannot be approximated by the lowest order eigenmode as 
traditionally done in SOI-MOSFETs. An empirical equation for the scaling length that depends 
on the geometry parameters is proposed. It is shown that even for a thick SiC>2 back oxide X e ff 
can be improved efficiently by thinner top oxide thickness, and to some extent, with high-k 
dielectrics. The model is then applied to self-consistent simulation of graphene nanoribbon 
(GNR) Schottky-barrier field-effect transistors (SB-FETs) at the ballistic limit. In the case of 
GNR SB-FETs, for large X e ff, the scaling is limited by the conventional electrostatic short channel 
effects (SCEs). On the other hand, for small X e f, the scaling is limited by direct source-to-drain 
tunneling. A subthreshold swing below lOOmV/dec is still possible with a sub-lOnm gate length 
in GNR SB-FETs. 
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I. INTRODUCTION 

Graphene has been attracting considerable attention as an alternative channel material to 
silicon in CMOS technology [l]-[6] because of its remarkable electrical properties [1] as well as 
due to the fact that an energy band gap can be induced under certain conditions in this normally 
gapless material [2]-[4]. Therefore, it may be possible to realize atomic-thin-body (ATB) 
transistors based on graphene. To date, GNR based transistor operation has been numerically 
simulated using the nonequilibrium Green's function formalism combined with 3D electrostatics 
[5]-[6]. A generalized scaling theory for the ATB field-effect transistor (ATB-FET) geometry 
however, has not been presented, which is different from the well established ultrathin-body 
(UTB) structures which have comparatively large body thicknesses with a well-defined 
semiconductor dielectric medium [7] -[9]. In this letter, a new generalized analytical solution is 
derived to calculate the electrostatic potential in the ATB-FET geometry which is applicable to 
transistors based on single-layer and few-layer graphene as well as conventional semiconductor 
materials with body thicknesses below ~ lnm [10]. On the other hand, carrier transport 
properties will depend on the individual bandstructure of the specific material. Here, we simulate 
GNR based Schottky-barrier FETs (SB-FETs) to explore their scalability dependence on 
geometrical parameters. 

II. APPROACH 

The structure of a double-gate (DG) ATB transistor is illustrated in Fig. 1(a). The atomic 
thin channel material with a length L g and a negligible thickness (assumed to be zero in the 
model) is placed in between the top and back oxides. The top oxide has a thickness of t tox and 
dielectric constant of e tox , and the back oxide is Si02 (£box = 3.9eo)with a thickness of hox- The 
potential boundary conditions are as shown in Fig. 1(a). 
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It should be noted that due to the negligible body thickness of ATB-FETs the parabolic 
approximation for the potential across the body cannot be used here as was done in original 
scaling studies [11]-[12]. Therefore, we use a procedure similar to that in [7]-[9], whereby the 
electrostatic potential of the ATB-FET can be written as y/{x, y) = v(x) + ul(x, y) + ur(x, y) , 
where v(x) is the long-channel solution that accounts for the top and bottom boundary conditions, 
as well as the self-consistent charge induced in the channel (therefore, the solution is valid in 
both subthreshold as well as above-threshold operation). ul(x, y) and ur(x, y) satisfy the source 
and drain boundary conditions, respectively, and can be expanded as: 



ul(x,y) = 



ur(x, y) = 



f sm[k n (x + t t J]smk n (L g -y)] 
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where the eigenvalues X n = l/k„ are determined by the implicit equation, 
s hox tan(k n t tox ) + s tox tsn(k n t box ) = . The coefficients Aj and An are calculated by integrating ul(x, 

0)*g n (x) and wipe, L g )*g n (x) from -t tox to tb ox , respectively, where g n {x) is constructed to be a set 
of corresponding conjugate functions to each term in ul(x, 0) that satisfies 

ttbox 

ul n (x,0)g m (x)dx = S n (note that g n {x) will simultaneously satisfy a similar condition with 

J— ttox 

each term of ur(x, L g ) as well). The calculated eigenvalues X n and coefficients A n l , A n r are shown 
in Fig. 1(b). Using the above series solution, an effective scaling length (X e ff) is extracted by an 
exponential fitting of the potential profile near the source region. Since the eigenmodes are 



4 



weighted by the expansion coefficients A„ and A n r , which do not necessarily peak at n = 1 (Fig. 
1(b)), l e ff cannot be approximated by the first eigenvalue (Ai) as has been traditionally done in 
SOI-MOSFETs [7]-[9]. For the geometry considered in Fig. 1(b) the first eigenvalue gives X\ ~ 
16nm, which shows a large discrepancy to the 2D series solution withA e ^ ~ 1.6nm. 

We use the above series solution to simulate GNR SB-FETs, considering a 2nm-wide 
GNR with an energy band gap of E g = 0.69eV [13], and source/drain SB height of <Pb = 0.1V. 
The Schottky barrier and direct source-to-drain (S->D) tunneling are both included on the same 
footing, whereby the WKB tunneling coefficient is calculated using the imaginary bandstructure 
of the GNR that treats the electron-hole nature of carriers realistically. Ballistic transport 
simulations are performed similar to [14] accounting for self-consistent electrostatics, while any 
edge-induced states in the band gap [15] and phonon scattering [16] are not considered. 

III. RESULTS AND DISCUSSION 
Figure 2(a) compares the extracted X e ff to the lowest order eigenvalue X\ as a function of t tox 
for different tbox values. It is clearly seen that X\ cannot correctly describe X e ff. Therefore, an 

empirical equation X = at l 

approximate X e ff, and this equation shows reasonable consistency for a large range of t tox and e tox 
values at thick SiC>2 box limit where a = 0.6 is the only fitting parameter (Fig. 2). It is seen in Fig. 
2(a) that X e ff is greatly improved by thinner t tox , and also modulated by tbox when it is below a 
certain thickness (< 50 nm). On the other hand, A^and X emp do not appreciably depend on ^ ox for 
thicknesses > 50nm. Figure 2(b) plots X e ff as a function of e tox for different t tox , and shows that k e g 
is also improved by high-k top oxide, but the improvement saturates at a certain point [5] [7]. 
More importantly, even for a thick Si02 back oxide, X e jf can be reduced below 2nm by thinning 
the top oxide thickness. 

5 



I _|_ £ box 



1 + 



box 
V £ tox J 



\Jbox J 



1 + 



' box 



is proposed to 



The transfer characteristics of GNR SB-FETs are shown in Fig. 3(a) for different gate 
lengths. For the calculated geometry, I on II ff> 10 is achieved at a supply voltage of 0.3V. 
Figures 3(b) and 3(c) show the band diagram and the energy-resolved current density in the off- 
state of the L g = lOnm and 15 nm devices, respectively. For this geometry which has X e g = 
1.45nm, when the channel length is L g > 15nm the subthreshold swing (SS) approaches the 
60mV/dec limit, but it degrades to 85mV/dec at L g = lOnm because of direct S->D tunneling 
(Fig. 3(b)). Here, we extract the SS values where logio(/z))-FGs is nearly linear in Fig. 3(a) at Vgs 
~ 0.05V. In order to further elucidate the scaling potential of GNR SB-FETs, Fig. 4(a) plots SS 
vs. L g , and Fig. 4(b) plots SS vs. L g l(nX e ff). Since the 2D potential at the middle of the channel 
near the thermionic barrier region depends on exp{-L g l22 e ff), previous works [7] [9] [12] have 
shown that SS degrades significantly for L g l(nX e ff) < 1.5 ~ 2 due to excessive thermal leakage. 

We, however, identify two distinct limiting conditions for SS degradation depending on the 
device geometrical parameters that are yet practically relevant: 1) for very small X e ff values, SS 
degradation is mainly due to direct S->D tunneling, while 2) at larger X e ff, SS degradation is due 
to conventional electrostatic short-channel-effects (SCEs) [7] [9] [12]. For thin t tox (lnm and 
2.15nm) in Fig. 4(a) with small X e g, SS degradation for L g below ~ 15nm is due to direct S^D 
tunneling [5], and not because of electrostatic SCEs. An SS < lOOmV/dec is still possible for sub- 
lOnm L g with small X e ff values, even though obtaining small off-state currents will be a challenge 
because of the direct tunneling limit. Nevertheless, narrower GNRs with larger band gaps may 
alleviate this problem [5]. For thick t tox (5nm and 8nm) in Fig. 4(b), on the other hand, with large 
Xeff, SS increases to ~100mV/dec when the conventional L g l(nX e ff) ~ 1.5 condition is met [7] [9] 
[12]. In other words, in the large X e ff limit, scaling is controlled by the electrostatic SCEs. 
Interestingly, when comparing Fig. 2 with SOI-MOSFET simulations in [9] for a similar device 
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geometry, the GNR channel has smaller X e g than the 5nm thick Si based UTB channel, showing 
the excellent scaling potential of GNR based ATB-FETs. 



IV. CONCLUSIONS 

A generalized 2D analytical electrostatic solution is developed for the ATB-transistor 
geometry. An empirical equation for the effective scaling length has been proposed which cannot 
be approximated by the lowest order eigenmode as traditionally done in SOI-MOSFETs. Even 
for a thick Si02 back oxide, the scaling length can be efficiently improved by thinner top oxide, 
and to some extent, by high-k dielectrics. It is seen that the scalability of GNR SB-FETs can be 
limited either by direct S->D tunneling, or conventional electrostatic SCEs, depending on the 
geometrical parameters. Compared to conventional UTB-FET geometries, excellent scaling 
potential of ATB-FETs based on graphene has been demonstrated. 
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FIGURE CAPTIONS 

Fig. 1 : (a) Schematic cross section of a double-gate (DG) ATB transistor. The channel material 
is buried between the top and the back oxides with zero thickness, (b) Eigenvalues X n 
and coefficients A n and A n r in the potential expansion (Eq. (1) and (2)) of a DG GNR 
SB-FET with t hox = 50 nm (Si0 2 ), t tox = 2.15 nm (high-k with e tox = 20e fl ), and V G s = V DS 
= 0.3 V. 

Fig. 2: (a) The extracted effective scaling length X e ff, the lowest order eigenvalue X\ and the 
empirical X emp as a function of t tox for different ttox- The back oxide is Si02 and the top 
oxide is high-k with e tox = 20eo. (b) X e ff and X emp as a function of e tox for different t tox . 
The back oxide is 50 nm Si02. 

Fig. 3: (a) Calculated transfer characteristics of DG GNR SB-FETs with different gate lengths 
at Kff = 1-45 nm. Band diagram and current density per energy in the off-state for a DG 
GNR SB-FET with (b) L g = 10 nm and (c) L g = 15 nm, showing that the off-state 
leakage is dominated by thermionic current at longer channel length (c), and increased 
dramatically by direct source-to-drain tunneling at short gate length (b) for this 
geometry (X e ff = 1.45 nm). 

Fig. 4: The subthreshold swing (SS) vs. (a) L g , and (b) L g IX e ff. For thin t tox (small X e fj), SS 
significantly increases for L g below 15nm due to direct source-to-drain tunneling, while 
for thick t tox (large X e fj), SS increases to lOOmV/dec when L g l(nX e ff ) ~ 1.5 due to 
electrostatic short-channel-effects (SCEs), as observed in conventional SOI-MOSFETs 
[7] [9] [12]. 
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Fig. 1. Q. Zhang, et al., IEEE Elect. Dev. Lett. 
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Fig. 2. Q. Zhang, et al., IEEE Elect. Dev. Lett. 




Fig. 3. Q. Zhang, et al., IEEE Elect. Dev. Lett. 
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Fig. 4. Q. Zhang, et al., IEEE Elect. Dev. Lett. 
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